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Topological stability of the edge states is investigated for non-Hermitian systems. We examine two 
classes of non-Hermitian Hamiltonians supporting real bulk eigenenergies in weak non-Hermiticity: 
SU{1, 1) and 50(3, 2) Hamiltonians. As an SU{1, 1) Hamiltonian, the tight-binding model on 
the honeycomb lattice with imaginary on-site potentials is examined. Edge states with Re_E = 
and their topological stability are discussed by the winding number and the index theorem, based 
on the pseudo-anti-Hermiticity of the system. As a higher symmetric generalization of SU{1, 1) 
Hamiltonians, we also consider SO{3, 2) models. We investigate non-Hermitian generalization of the 
Luttinger Hamiltonian on the square lattice, and that of the Kane-Mele model on the honeycomb 
lattice, respectively. Using the generalized Kramers theorem for the time-reversal operator Q with 
0^ = -1-1 [M. Sato et al., arXiv:1106.1806 , we introduce a time-reversal invariant Chern number 
from which topological stability of gap less edge modes is argued. 
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I. INTRODUCTION 

Recent developments of topological insulators have led to much interest in topological phases of matter. Murakami, 
Nagaosa, and Zhang predicted spin-Hall effect [l|,[4| by re-examining Haniiltonians describing hole-doped semiconduc- 
tors with spin-orbit coupling presented by Luttinger (the Luttinger Hamiltonian) 3] . After the theoretical predictions 
and experimental observations [3, Q, the spin-Hall effects were generalized to insulating systems [y|- Subsequently, 
the quantum version of the spin Hall effect, the quantum spin Hall effect, was introduced by Kane and Mele as a 
model with time-reversal symmetry in graphene [Tj and independently by Bernevig and Zhang as two-dimensional 
semiconductor systems with a uniform strain gradient [8| . Remarkably, quantum spin- Hall effects were experimentally 
observed in 2D CdTe/HgTe/CdTe quantum well by Konig et a/.^i following the theoretical predictions by Bernevig, 
Hughes, and Zhang |l3l. 

The appearance of the topologically protected gapless edge states within the bulk gap is a manifestation of the 
topological insulator. The number of such gapless edge modes is specified by topological invariants. In Ref . [ll| . 
gapless edge states carrying spin currents were found when the spin is conserved. They are characterized by spin 
Chern numbers introduced by Sheng et a/.[l2|. The spin Chern number is an extension of the Chern number for 
quantized Hall conductivity [13|, Il4| to quantized spin Hall conductivity in time-reversal symmetric systems with spin 
conservation. Kane and Mele introduced Z2 invariants that distinguish topologically non-trivial phases from trivial 
ones in time- reversal symmetric systems [TJ, |l5| even when the spin is not conserved. The above-mentioned gapless 
edge states in quantum spin Hall insulator, dubbed as the helical edge modes [l3|, are known as a consequence of the 
Kramers theorem of the time-reversal symmetry Q with Q^ = —1. In contrast, time-reversal symmetric systems with 
Q^ = -1-1 have not been investigated so far in the context of topological insulators, since they are considered to be 
irrelevant to topologically protected gapless edge states. 

Recently, the present authors have shown that the (generalized) Kramers theorem follows even for 9^ = +1 in a 
class of non-Hermitian systems (provided the metric operator is anticommutative with the time- reversal operator) [17]. 
Therefore, a natural question arises: Do topologically protected edge modes also appear in such non-Hermitian models? 
The main purpose of the present work is to give an answer to the question. We demonstrate numerical calculations 
and provide topological arguments for the stability of edge modes in non-Hermitian Hamiltonians. In particular, we 
inve stig ate lattice versions of the 5*^7(1, 1) and SO{3, 2) Hamiltonians studied in Ref. [l3|- (See also the related works, 
[la . [l9| in higher dimensional quantum Hall effect, and 20l-|22| in PT symmetric quantum mechanics.) As a lattice 
realization of 5J7(1, 1) model, we consider the tight-binding model on the honeycomb lattice with imaginary on-site 
potentials. For 50(3, 2) model, we investigate non-Hermitian generalization of the Luttinger Hamiltonian [y, [ll| on 
the square lattice. We also argue a non-Hermitian generalization of the Kane- Mele model[7|, [l5|, where the hopping 
integrals are asymmetric due to non-Hermiticity. 

Non-Hermitian systems play important roles in physics. For instance, a non-Hermitian system with disorder, known 
as the Hatano-Nelson model, has been studied in the context of localization-delocalization transition in one and two 
dimensions 23 25]. The model simulates the depinning of flux lines in type-II superconductors subject to a transverse 
magnetic field [2^. We also briefiy discuss possible relevance to the model. 

The paper is organized as follows. In SecHH we introduce the SU{1, 1) and 50(3, 2) Hamiltonians and their basic 
structures. In Sec lIIIl a lattice version of the 5f/(l, 1) model on graphene is explored, and the stability of edge states 
is discussed on the basis of topological arguments. We further investigate the 50(3, 2) model on the square lattice 
in Sec lIVI A topological number for the non-Hermitian system is introduced in order to account for the stability 
of edge modes under small non-Hermitian perturbation. In Sec|Vl a non-Hermitian version of the Kane-Mele model 
is presented and the topological property of the model is also discussed. Section IVII is devoted to summary and 
discussions. 

II. SU{1, 1) AND SO{3, 2) HAMILTONIANS 

In this paper, we mainly consider SU{1, 1) and 50(3, 2) models, which are non-Hermitian generalizations of SU{2) 
and 50(5) models, respectively. Here, we briefly discuss their structures with emphasis on the relations to split- 
quaternions along Ref. Tr]. 

It is well known that the quaternion algebra is realized as 2 x 2 unit matrix, I2, and (i times) Pauli matrices, iaa 
{a = 1, 2, 3). The SU{2) Hamiltonian is constructed by the Pauli matrices aa, 

H = diai + d2a2 + dsas, (1) 

where da (a = 1, 2, 3) are real. The eigenenergies are. 



E±=±Jdl + dl + dl. (2) 



Similarly, the split quaternion algebra is realized as the 2x2 unit matrix, I2, and (i times) SU{1, 1) Pauli matrices, 
(— (Ti, — (72, iCTa), and the SU{1, 1) Hamiltonian is constructed by such SU{1, 1) Pauli matrices: 

H = di<Ti + d2(J2 + id^as. (3) 

The eigenenergies are. 



E±=±^dl + dl-dl (4) 

The 5*0(5) gamma matrices, Ta (a = 1, 2, ... , 5), are defined so as to satisfy the Clifford algebra {Fa, Ff,} = 26ab- 
They have the quaternionic structure since their off diagonal components are given by the "quaternions" 

^. ^ (-L. 7) ■ r» - {X '7) ■ -3 - i-l 7) . ^' - C ;?) . ^= - {]' X) - (^) 

With the 5*0(5) gamma matrices, the 50(5) Hamiltonian is given by 

F = diFi+d2F2+d3r3 + d4F4 + d5r5, (6) 

where da (a = 1, 2, . . . , 5) are real. The eigenenergies are, 



E±=±^dl+dl + dl + di + dl (7) 

each of which is doubly degenerate. Such double degeneracy is understood as a consequence of the Kramers theorem 
since the 5*0(5) Hamiltonian is invariant under time-reversal operation 0^ = — 1 with 

where K denotes the complex conjugation operator. Similarly, the 5*0(3, 2) gamma matrices are introduced as the 
matrices whose off-diagonal blocks are given by the split-quaternions. Hence, the 5*0(3, 2) gamma matrices are given 

by 

-ai\ fO ~a2\ ( ia:A fO IsA fh \ 

which are iFi, iF2,F3,F4,F5 with the 50(5) gamma matrices ([5]). From the 50(3,2) gamma matrices, we construct 
the 5*0(3, 2) Hamiltonian as 

H = idiFi + id2T2 + dsFa + d4F4 + d^T^. (fO) 

The eigenenergies are. 



E±=±^^d'3 + di + di-di-di, (ff) 

each of which is doubly degenerate. The 5*0(3, 2) Hamiltonian is invariant under time-reversal operation 0^ = +1 
with 

As discussed in Ref . \ll\ . such double degeneracy is a consequence of the generalized Kramers theorem for 6^ = +1. 
The metric operator that satisfies rjHrj^^ = W is given by 

Apparently, rj is anticommutative with 0+: {77,0+} = 0. 
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FIG. 1: The honeycomb lattice. 




FIG. 2: The honeycomb lattice with zigzag edges at m = 1 and m, = 8 [Ly — 8) along the i-direction. 



III. 5(7(1, 1) MODEL AND WINDING NUMBERS 



In this section, we investigate topological stability of the edge states of the SU{1, 1) model. As an example, the 
tight-binding model on the honeycomb lattice with imaginary on-site potentials is considered. It will be found that 
edge states with KeE — are robust under weak non-Hermiticity. We will also argue that their topological stability 
is guaranteed by topological reasoning. 



A. Graphene with imaginary sublattice potential 

As a model with the SU{1, 1) Hamihonian, we consider the following tight-binding model on the honeycomb lattice: 

ff-t^(clc, +H.c.)-i^A,c|c„ (A, >0). (14) 

The first term is the nearest-neighbor hopping. The second term is imaginary on-site sublattice potentials which make 
the system non-Hermitian. Here, A^ = Xa and A^ for closed and open circles in FiglU respectively. The Hamiltonian 
(ITi)) is rewritten as 

H = tJ2icU+^■C.)+^XvY.^^cU-^^^^Y.4c^, (15) 

where Ay = ^''^'^•^ and ^i = -1-1 (—1) for closed (open) circles in FiglTJ Since the — z '^"^'^'' ^^ c^^Ci term only shifts 
the origin of the energy, we neglect it in the following: 

H = tY.^c\cj + H.c.) + iAy ^e.clc,. (16) 

In the absence of edges, (J16p reduces to the following SU{\, 1) Hamiltonian in the momentum space: 

iXv 2icos(%^) -l-te*''"' 



H{k)=\ 7 .^,„„ —y ^ >''- - , (17) 

V 2t cos (^) + te-*^^ -~iXv ) 

which is obtained by the Fourier transformations, 

k 

c.+ci = 0(„+i.™+i)-Ee'*'^"^'"^''''™+'Vfc- (18) 

k 

Here, V'fc and (pk are the annihilation operators in the momentum space corresponding to the closed and the open 
circles in Figll] respectively, kx (ky) the momentum in the x (y) direction, a the lattice constant of the honeycomb 
lattice, and d = ay/\/3. The bulk spectra of the system are obtained by diagonalizing the Hamiltonian (|17l) . 



E±ik) = ±^ 



, k^-aX .'/3ky 
t-^2^cos ( ^ I e*^ 



X'v (19) 



Let us now examine the edge state in this system. If we make zigzag edges along the x direction as shown in FiglJJ 
a zero energy edge band appears in addition to the bulk bands (|19p : For the Hermitian case Ay = 0, it has been 
known that the edge state with E — appears for 27r/3 < akx < 47r/3 as shown in Fig|3] [23 • We find here that the 
zero energy edge state persists even in the presence of weak non-Hermiticity Ay. (See Figs|4]and[5l) We illustrate 
the real part and the imaginary part of the energy bands as functions of akx in Fig|4] and FiglS] respectively. It is 
clearly seen that the edge state with Re -E = persists in the region 2tt/3 < ak^ < 47r/3 for small Ay. 

B. Winding number and generalized Index theorem for non- hermitian system 

From the bulk-edge correspondence, the zero energy edge state suggests the existence of a topological number 
responsible for the edge mode. In this section, we will see that this is indeed the case and the edge state is characterized 
by the one-dimensional winding number. 

1. Basic property of edge states 

In order to make our arguments concrete, we consider the semi-infinite SU{1, 1) model on m > 0. We perform the 
inverse Fourier transformation of H{k) in (1171) with respect to ky and denote the resultant Hamiltonian as H{kx)m,m' ■ 
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FIG. 3: The energy bands of the SU{1, 1) model p6|) with zigzag edges along the a;-direction for t = 1.0 and \v = (Hermitian 
case). Here a is the lattice constant and k^ the momentum in the s-direction. A zero energy edge state with fiat band appears 
for 27r/3 < ak^ < Atv/3. 




FIG. 4: The real part of the energy bands of the SU{1, 1) model (|16p with zigzag edges along the a;-direction. We show the 
results for t = 1.0 and various values of non-Hermitian parameter Xv- Here a is the lattice constant and k^ the momentum in 
the a;-direction. For weak non-Hermiticity [(a) and (b)], the edge state with KeE — survives around ak^ ~ n. For large Xv 
[(c)-(f)], no edge state with KeE = appears. 



The energy spectrum of the system is given by 



(20) 



with 



\u{k^,m)) = 0, 



(21) 



for TO < 0. In the following, we consider k^ as a parameter of the system and treat the model as a one-dimensional 
system along the y-direction. We also restrict our argument on k^. where the bulk energy gap is open in the real part 
of the energy. The zero energy edge state satisfies ^U\i and (I^T]) with KeE(kx) = 0. 
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FIG. 5: The imaginary part of the energy bands of the 5)7(1, 1) model (|16p with zigzag edges along the s-direction. We plot the 
results for t = 1.0 and various values of the non-Hermitian parameter \v- Here a is the lattice constant and k^ the momentum 
in the a;-direction. For weak non-Hermiticity [(a) and (b)], there is no imaginary part in the bulk except near the gap closing 
points at ak^ = 27r/3 and akx ~ 47r/3. The edge states with KeE — [Fig|4] (a) and (b)] support the imaginary part in their 
energy. For large Av [(c)-(f)], the bulk states support the imaginary part. 
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FIG. 6: (a) Topologically protected state with ReS = 0. (b) Topologically trivial states with Re£ = 0. Here, the closed (open) 
circles at Rei? = represent states satisfying a-i\uf) = + \u) (|27p {a-i\uf) = — \u) (|28p '). In the latter case [(b)], the state can be 
non-zero mode since n\ — n^L — 0. 



As was shown in Ref.[l7|, the basic symmetry of the SU{1, 1) Hamiltonian is pseudo-anti-Hermiticity given by 

a3H{k)^a-s = -H{k), (22) 

which also implies the pseudo-anti-Hermiticity of the Fourier transformed one, 

0-3[H{kxy]m.Tn'Cr3 = - H{k^)ra,7n' ■ (23) 

The pseudo-anti-Hermiticity is a key ingredient of our topological argument. 



Let us denote the right eigenvectors and the left eigenvectors oi H{kx)m,m' as \u{kx, m)) and \u{kx, rn))), respectively: 

^ H{kx)m,rn' ["(^x, m')) = E{kx) \u{kx,m)) , 

Y,[H{kxyUra'\u{kx,m'))) ^ E{kxy\u{kx,m))), (24) 

m' 

where |u(fca;,'7i)) and |w(fcj:,TO))) are normalized as 

^{{u{kx,m)\u{kx,Tn)) = '^{u{ka;,m)\u{kx,Tn))) = 1. (25) 



From the pseudo-anti-Hermiticity (j23l) . we find that 

^H{kx),n,m'0-3\u{kx,m'))) = -E{kx)*(J3\u{kx,m))). (26) 

Therefore, the eigenstate |u(fea;,'7i)) with eigenenergy E{kx) always comes in pairs with the eigenstate cr3|M(A:2;,m))) 
with eigenenergy —E{kx)*- When 'R,eE{kx) ^ 0, |u(fca;,m)) and (T3|u(fca;,m))) are independent of each other, since 
they have different energies. On the other hand, for edge states with I{BE{kx) = 0, they are not always independent. 
Actually, by choosing a proper basis, they can be related to each other as 

az\u{kx,m))) = +\u{kx,m)) , (27) 

or 

(Jz\u{kx,m))) ^ -\u{kx,m)) , (28) 

where the overall phase factors of the right-hand sides of (l?f|) and (1^51) are restricted to ±1 due to the normalization 
condition (p5)) . We denote the number of the edge states with Re£'(fcx) = satisfying (l27t and (|28)) as n*]. and n*L, 
respectively. 

Here we can show an important property of the edge states: The difference between n*]_ and nP_ does not change its 
value against perturbation preserving the pseudo-anti-Hermiticity. The reason why rP^ — nP_ is conserved is as follows. 
For IieE{kx) ^ 0, \u{kx,m)) and a3\u{kx,m))) are independent, so we can construct another basis from them as 

\u±{kx,m)) = —={\u{kx,m))) ± a^ \u{kx,m))), 

\u±{kx,m))) = —={\u{kx,m))±(T3\u{kx,m)))), (29) 

with 

{{u±{kx,m)\u±{kx,m)) ^ 1, {{u±{kx,m)\u^{kx,m)) ^ Q. (30) 

Since these states \uj^{kxTm)) and \u^{kx,m)) satisfy (|27l) and (1^51) . respectively 

(T3\u±{kx,m))) ^±\u±{kx,m)) , (31) 

we have the same number of states with signs -I- and — in the right-hand side of (I5T|) for IleE{kx) ^ 0. This means 
that vP^ — n° can not change adiabatically: By small perturbation, some of the edge states may acquire a non-zero 
real part of the energy. If this happens, however, they must be paired with opposite sign in the right hand side of 
(PIjl . As a result, the difference between n^ and nP_ does not change at all. In FiglHl we illustrate two different edge 
states with Re£' = 0, one is topologically protected from the above argument and the other topologically trivial. 

2. Winding number 

Now we introduce a bulk topological number relevant to the present edge mode. To do this, consider the eigenequa- 
tion for H{k) in the momentum space: 

H{k) |u„(fe)) = En{k) \ujk)) , i7(fc)t|u„(fe))) = En{k)*\u^{k))), (32) 



where n denotes the index labehng different bands. We assume that the system is half-filhng and that we have a gap 
in the real part of energy around KeE{k) = for a fixed value of kj.. 

From the pseudo-anti-Hermiticity (|22|) . one can say that if |u„(fe)) is a right eigenstate of H{k) with Re-E„(fc) > 0, 
(T3|u„(fc)) is a left eigenstate of H{k) with Rei?„(fc) < 0. In the following, we use a positive (negative) n for the state 
with ReEn{k) > (Rei?„(fc) < 0), and set the relation 

|u_„(fc)» = as \unik)) , E^„{ky = -Enik). (33) 

The topological number is constructed from the non-Hermitian generalization of the projection operators Vi{k) 
and V2{k) for occupied bands J27l.[28l|: 

Pi(fe) - Y. i"»(^))) ("»wi ' ^2(fe) = Y. i""(^)) (K(fc)i, 

(Vi{k)y = V,{k), {V2{k)f = Mk). (34) 

From 7^1 (fc) and 'P2{k), we define the following Q matrix: 

Q{k) = 1- [Vi{k)+V2{k) 

\n>0 n<0 

+ E l"n(fe)) (K(fc)l - E l^"(^)) (("«(fe)l ) ' (35) 

n>0 n<0 / 

where the completeness relation J^n Wn{k))) (u„(fc)| — J^n Wn{k)) ((u„(A;)| = 1 was used. One can immediately show 
that the matrix Q{k) is Hermitian: 

Q(fc)^ = Q(fc). (36) 

From (p3| and ((35)) . we find that Q(fc) and 0-3 anticommute: 

{Q(fc),a3} = 0. (37) 

Since 0-3 is diagonal, ([37| implies that Q(A;) is off-diagonal and can be expressed by a complex function q{k) as 

The topological number relevant to our model is the one-dimensional winding number wid defined as 

wu{k.) ^ T^ f q{k)-^^q{k)dky. (39) 

By using polar coordinates qik) — |g(A;)|e'"'-'°\ we obtain, 

1 /'27r cj /'jU~\ 1 /'27T cj 

= iV, (TV : integer), (40) 

where we have used the periodicity of q{k) with respect to ky to obtain a quantized value A''. The bulk-edge corre- 
spondence implies that if the winding number wid{kx) is nonzero, there is a zero energy state (in the real part of the 
energy) on the boundary. We will confirm this numerically in Sec lIII Cl 
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3. Generalized Index theorem 

In Sec lIIIB"2l we argue that the non-zero wid{kx) implies the existence of edge states with Re£' — 0. At the same 
time, in Sec lIIIBTl we find that the non-zero nP^_ — n^ ensures the robustness of the existence of edge states with 
Re-E = 0. Therefore, it is natural to identify these two quantities Wid(fcx) and n^_ — nP_. Since there is a sign ambiguity 
for the identification, two possible relations are suggested in the form of the index theorem, 

n"_-nl= wid(fcx), (41) 

or 

nl-n°_= wid(fcx). (42) 

Here note that there are two possible choices of the edge of the system, i.e., the surface of the semi- infinite system on 
y > or that on y < 0. The two possible choices of the surface correspond to two possible equalities (|4ip and (j35]). 

In this paper, we will not prove the generalized index theorem (|4T|) and (|42|) . In Ref.pO!], one of the present authors 
proved a similar generalized index theorem for zero energy edge states in systems with chiral symmetry. We expect 
that the generalized index theorem in this case can be proved in a similar manner. 

C. Application to SU{1, 1) Hamiltonian 

Here we present explicit calculations of the winding number wi(j(fcj;) for our SU{1, 1) Hamiltonian (|17p . Since we 
treat k^ as a fixed parameter of the system, it is convenient to write (1171) as 



with the identification 



Let us write (1551) as 



^ = Xy, v^2tcos('^), v'^t, k^y^. (44) 



Q(fc) = Q'(fc) + Q'{k)\ (45) 



where 



\n>0 

Then, straightforward diagonalization of (|33|) gives 



2 , 



Q^w- . . '.. A -2-..''^''^ 



1 



2^|f -hw'e»'=|2-72 
From (SSI) and (HTl), we obtain 



i?(fe)^ (47) 



Q(fc) = ] =[H{k)^+H{k)] 



1 / v + v'e 



y/\v + v'e'''\'^-j^ \v + v'e '^ 



I „ik 



(48) 



Therefore, the winding number widikx) is evaluated by ([39| with qik) = , . [v + v'e ). Since only the 

{v + v' e^^) part contributes to the winding number, we have 

wuikx) = — y z^^j-^dk. (49) 



11 

Furthermore, changing variable as z = e*'^ leads to 

v' f 1 

wid{kx) = 7-- (p — ; — —dz, (50) 

ZTTl J V + V'Z 

where the integral is taken over a unit circle \z\ — 1. Using the residue theorem, we reach the final result, 

widik,) = 1, i\v'\ > \v\), 

= 0, {\v'\ < \v\). (51) 

Therefore, the bulk-edge correspondence predicts the existence of the edge state with KeE = if \v'\ > \v\. 
In terms of the original parameters of the Hamiltonian ([T7|) . the inequality \v'\ > \v\ means 



cos ; ^2" 



< 1. (52) 



Thus the edge state with Re£' = is predicted in a region of 2tt/3 < ak^ < 47r/3, provided k^ also satisfies 
(2t\ cos (%^) I ^ ^) > A^ so that the bulk energy gap is open. This prediction is clearly confirmed in Figs. mandO 

IV. 50(3, 2) MODEL AND TIME-REVERSAL INVARIANT CHERN NUMBER 

In this section, we examine SO{3, 2) Hamiltonians as a higher symmetric generalization of the SU{1, 1) ones exam- 
ined in Sec. IIIIl As a concrete example, we consider a non-Hermitian generalization of the Luttinger Hamiltonian |6l. 1 1 1| 
on the square lattice. We find that gapless edge states exist under weak non-Hermiticity. To explain the topolog- 
ical origin of the gapless edge states, we introduce a time-reversal invariant Chern number inherent to a class of 
non-Hermitian Hamiltonians, based on discrete symmetry of the system. From the bulk-edge correspondence, the 
time-reversal invariant Chern number ensures the existence of the gapless edge states in the non-Hcrmitian system. 

A. General Hamiltonian with time-reversal symmetry and pseudo-Hermiticity 

We first derive a general 4x4 non-Hermitian Hamiltonian with pseudo-Hermiticity |30l . |3l| , which is invariant under 
time-reversal symmetry 0+ with Q^ = +1 as well. Such a Hamiltonian belongs to one of the 43 classes of random 
matrix classification for non-Hermitian systems presented in Ref . p33| . Gapless edge states obtained in the following 
are expected to be robust against disorder in the same class. 

A general 4x4 complex Hamiltonian H{k) can be represented by a linear combination of the identity matrix, 5 
gamma matrices F^, and 10 comnuitators Tab = [^a, rb]/(2i) (a, 6 = 1, 2, • • • ,5): 

5 5 

H{k) = ho{k)+J2ha{k)ra+ Y. hab{k)Tab, (53) 

a=l a<6=l 

where ho{k), ha{kys and hab{kys are complex functions of k. We adopt the following representation of the gamma 
matrices: 

■ I2 A T. f ^ -«l2 \ ^ f fT"-2 



^^-[i\ oj' r.= ^i j' '^^ = [ -."-j' (54) 

where a = 3, 4, 5 and a^ {p = 1, 2, 3) are the Pauli matrices [39{ . 

Now we consider time-reversal symmetry. The time-reversal operator Q is represented as 

Q^U -K, (55) 

where C/ is a unitary operator, and K is the complex-conjugate operator. The square of a time-reversal operator is 
either +1 or —1; 

e^=±l. (56) 

In the following, we focus on the time-reversal symmetry 9+ with O'^ = +1. For later use, it is convenient to choose 
U in 0+ as [/ = riF4. Imposing the time-reversal invariance on H[k), 

e+iJ(-fe)e7^ = i/(fe), Q+^TiVi-K. (57) 
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we obtain 

/lo(fc) = ho{-k)*, /l(i,2)(fe) = -/l(l,2)(-fe)*, /l(3,4,5)(fe) = ^(3,4,5)(-fc)*, 

^(12,34,35,45) C*^) = ^''■(12,34,35,45) (^''^) i ^(13,14,15,23,24,25)1''') = ^1(13,14,15,23,24,25) ("''') • (58) 

In addition to the time- reversal invariance, we impose pseudo-Hermiticity on H{k), 

r]H{k)^7^-^ =H{k), v = i^2Tu (59) 

with Hermitian matrix "q. rj is cahed the metric operator. Here we have supposed that the metric operator rj anti- 
commutes with the time-reversal operation, 

{e+,77}=0, (60) 

which restricts the allowed form of rj. From a proper unitary transformation and rescaling, we can always take the 
basis in which rj is given by 77 = ir2ri. Equation ([59]) leads to 

ho{k) = ho{k)*, /l(l,2)(fc) = -/i(l,2)(fc)*, /l(3,4,5)(fc) = ^(3,4,5)(fc)*, 

''•(12,34,35,45) C*!) = ^(12,34,35,45) (fe)* j ^(13,14,15,23,24,25) (fe) — — ^(13,14,15,23,24,25) C*!)*- (61) 

Combining ((58| and ((6T|) . H{k) is written as 

F(fc)=ao(fc)+^ ^ 6^(fc)r^+ Yl «p(fc)r/x 

Ai=l,2 p=3,4,5 

+ ^ a^,(fe)r^, +^ Yl b^,{k)T^,, (62) 

^1^=12,34,35,45 a'!^=13, 14, 15,23,24,25 

where ao{k), ap(A;)'s, 6^(A;)'s, ap^(fc)'s, and h^^{kys are real functions of fc, and satisfy 
ao(fc) == ao(-A;), 6(i,2)(fc) = fc(i,2)(-fc), a(3,4,5)(fc) = a(3,4,5)(-fc), 

^(12, 34,35,45) e*!) = ~a(12, 34, 35,45) ( — ''5)7 ^(13,14,15,23,24,25) (fc) == ^^(13,14,15,23,24,25) ("fc)- (63) 

B. 50(3, 2) Luttinger model and edge state 

In order to have real eigenenergies in non-Hermitian systems, PT symmetry plays crucial roles as pointed out by 
Bender et a/. |33l-l35l| . Following the arguments by Bender et ai, we impose the inversion symmetry on H{k) in (|62|) : 

H{-k)=H{k). (64) 

From (|62]), dBS]), and (jMl), we obtain the following 5*0(3,2) Hamihonian [l7 J: 

ff(fc)=ao(fc)+^ ^ 6^(fc)r^+ Y. aA^e^)^/^' (65) 

Ai=l,2 p=3,4,5 

with real even functions ao(A;), ap(fc)'s, and fep(A;)'s. The eigenenergies are 



£;± (fc) = ao(fc) ± Va3(fc)2 + a4(fc)2 + 05(^)2 - 6i(fc)2 - &2(fc)2, (66) 

each of which is doubly degenerate. As is expected from PT symmetry, these eigenenergies are real if a-i{k)'^ + ai{ky' + 
a^{kf>bi{kf+b2{kf. 

To realize such an 5*0(3, 2) model, we generalize the 5*0(5) Luttinger Hamiltonianjy, [ll| on the square lattice into 
a non-Hermitian form: 

H{k) = e{k) + V Y daik)ra + iV Y '^-^-^ (67) 

a=3,4,5 a=l,2 

where 



e(fc) ~ i(— 2 cosfcj; — 2 cos ky + Cg), c?3(fc) = — vSsinfcj: sinfcy, 

d4(A;) = \/3(cos/c2; — cosfcj^), d5(A;) = 2 — e^ — cosfc^^ — cosfcy, (68) 
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di-d2-0 




Ky 

FIG. 7: The energy bands of the 50(3, 2) Luttinger Haniiltonian (|67|) with open boundary condition in the x-direction for 
t — 1.0, V — 4.0, es — 0.5, and di = d2 = (Hermitian case). Here ky is the momentum in the y-direction. Two gapless helical 
edge states appear on each edge. 



with k = (kx, ky). The imaginary unit i in the last term in the right hand side of (|67p is absent in the original 5*0(5) 
Luttinger Hamiltonian. For simplicity, we take di and ^2 to be real constants. Here, Cg = (fc^) where ( ) represents 
the expectation value in the lowest band [401. The eigenenergies are 



E±{k) = e(fe) ± VJdsikf + diiky + d^ikf -dl- dl 



(69) 



each of which is doubly degenerate. 

Now we examine edge states of this model. We put the system on a cylinder with the periodic boundary condition 
in the y direction and open boundary condition in the x direction, and study the quasiparticle spectrum numerically. 

First, we illustrate the quasiparticle spectrum for the Hermitian case with di — d2 — 0. The energy bands in this 
case are shown in Fig. [T] In the gap of the bulk bands, there exist doubly degenerate edge bands. The existence of 
the gapless edge bands is explained by so-called the spin Chern number [Til. When di = d2 = 0, in addition to the 
Q^ = +1 time-reversal symmetry ([57|). the Hamiltonian (|67|) has the following Q^ — —1 time-reversal symmetry: 

(70) 



e-H{-k)QZ^ =H{k), Q^=T2Ti- K. 
Moreover, since H{k) — H^{k) for di = d2 = 0, the pseudo-Hermiticity ([5^ reduces to 

[iJ(fc),7?]=0. 



(71) 



Thus in this special case, ry becomes a conserved quantity. Indeed, the operator r/ is identified with the pseudo-spin 
operator Sz in the Luttinger model [ll| whose eigenvalue is either +1 or —1. In Ref . |ll| . the spin Chern number 
Gs was defined by using the pseudo-spin Sz and the time-reversal symmetry 9_, and it was shown that Cs = 2 if 
< e^ < 4. Therefore, the existence of the gapless edge bands is ensured by the spin-Chern number. 

Let us now examine the non-Hermitian case. We show the real part and the imaginary one of the energy bands 
as functions of ky for various ((ii,(i2)'s in Figs. [S] and El respectively. Interestingly, the gapless edge bands persist 
in the real part even in the presence of the non-Hermiticity {di,d2) [Figs 15] (a) and (b)], although they have a small 
imaginary part at the same time [Figsl9l (a) and (b)]. Since the non-Hermiticity breaks both of the time-reversal 
invariance for Q- and the rotational invariance for Sz, such gapless modes cannot be explained by the spin Chern 
number. The topological origin of the gapless edge mode will be discussed in the next section. If we further increase 
(di,d2), the bulk gap closes, and only remnants of the edge states are observed [Figsl8l(c)-(f)]. 



Time-reversal invariant Chern number 



In the previous section, we showed that there exist gapless edge states in the SO{3, 2) non-Hermitian model, 
this section, we will discuss the topological origin of these gapless edge states. 



In 
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di^om di=om (b) 



rfl=0.1 d2=01 




FIG. 8: The real part of the energy bands of the 50(3, 2) Luttinger Haniihonian (|67|l with the open boundary condition in the 
x-direction. We plot the resuhs for t = 1.0, V = 4.0, Ea — 0.5, and various values of the non-Hermitian parameters {dijd^). 
Here ky is the momentum in the y-direction. For weak non-Hermiticity [(a) and (b)], two gapless hehcal edge modes survive. 
When the bulk gap in the real part of the energy closes [(c)-(d)], the topological number Ctri is not well defined. In this 
parameter region, only remnants of the edge states are observed. 



1. Generalized Kramers theorem 

First of all, we examine general properties of eigenstates for non-Hermitian Hamiltonian. Let us denote the right 
eigenvectors and the left eigenvectors of a non-Hermitian Hamiltonian H{k) as \un{k)) and |u„(fc))), respectively: 



H{k) |u„(fc)) = E,,{k) \u^{k)) , iJ(fc)t|u„(fc))) = -B„(fc)*|w„(fc))), 
where |u„(fc)) and |m„(A;))) are normalized as 

{{Um{k)\Un{k)) = (u™(fc)|u„(fc))) =(5m„. 

The eigenstates |w„(fe)) and |M„(fc))) satisfying ((73)) are called as the bi-orthonormal basis. 
Now suppose that H{k) is pseudo-Hermitian: 

i7(fc)^ ^r]H{k)r]-\ 

From the pseudo-Hermiticity, the first equation of ([7^ is rewritten as 

H{k)^7j\Un{k)) ^ Er.{k)7j\Ur.{k)). 

Therefore, ri\un{k)) is an eigenstate of H{k)^ and it can be expanded as 

'q\Un{k)) = ^ \Ur,^{k)))Cr,^n{k)■ 

rn 

Applying {um{k)\ from the left, we obtain 



(72) 
(73) 

(74) 
(75) 
(76) 



Cmn(fc) = {u,n{k)\r]\Un{k)) . 



(77) 
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FIG. 9: The imaginary part of the energy bands of the SO{3, 2) Luttinger Hamiltonian H67p with the open boundary condition 
in the s-direction. We plot the results for t = 1.0, V = 4.0, Cs — 0.5, and various values of the non-Hermitian parameters 
(di,d2). Here ky is the momentum in the y-direction. For weak non-Hermiticity [(a) and (b)], there is no imaginary part in 
the bulk bands. Only the edge states with KeE = [Fig[S] (a) and (b)] support the imaginary part in their energy. When the 
bulk gap in the real part of the energy closes [Fig[8l^c)-(d)], the bulk states also support the imaginary part as well. 



Since rj is Hcrmitian, Cm„(fc) is also Hermitian with respect to the index m and n. Thus it can be diagonalized by a 
unitary matrix G: 



'^Gmi{ky cik{k)Gkn{k) = A„i(fc)(5„i„, 



(78) 



Ik 



with real A„(fe). The eigenvalue A„(fc) is not zero since Cmn{k) is invertible. Taking the following new bi-orthonormal 
basis 



|0„(fc)) = Y. \Um{k))Gmn{k)/^J\\n{k)l |0„(fc))) = ^ \u„,{k)))G,nn{k)^J\\Jk)l 



(79) 



with {(j)ni{k)\(j)n{k))) = (((/)m(fc)|</)„(fc)) = S^n, we havB 

77|</)„(fe))=sgn(A,„(fe))|</)„(fc))). 



(80) 



Because of the continuity of the wave function, sgn(A„(fc)) does not change in the whole region of the momentum 
space. Thus, the states in the new basis are classified into two, i.e. states with 



and those with 



v\Mk)) = |0„(fc))), 



r,|0„(fc))=-|</>„(fe))). 



(81) 



(82) 



Here we should note that the new bases |^„(fc)) are no longer eigenstates of H{k) unless -E„(fc) is real. Indeed, the 
right hand side of ([7^ mixes the eigenstates with En{k) and those with En{k)* . However, the mixed states have a 
common real part of the energy. Therefore, even if the bulk energy Enik) is not real, the structure of the real part 
of the eigenenergy remains the same in this basis. 
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If H{k) is invariant under the time-reversal symmetry 0+ with 9^ = +1 and {Q+,r]} = 0, we have an additional 
structure, i.e. the generalized Kramers degeneracy [i7|: Let |(/)„ (k)) be a wavefunction satisfying (|8ip . 

77|0W(fc)) = |4+)(fc))>, (83) 

with 

{<j>i^HkMl+\k))) = ((0(+)(fc)|0(+)(fc)) = 5„„. (84) 

The generalized Kramers partner |(^„ (fc)) is given by 

|0i-)(fc))=e'«"(^)e,?-Vi+n-fe))), (85) 

with a phase factor 9n{k), and the corresponding left state is 

|0i-)(fc)))=e^«"('=)e7?-i|^W(~fe)), (86) 

with 

(0(,7)(fc)|0(-)(fc))) = ((0(,;)(fc)|0(-)(fc)) = <5„„. (87) 

The generalized Kramers pairs \4>n (k)) are independent of each other [l3|- Actually, it is found easily that at time- 
reversal invariant momenta k = Ti satisfying ~Ti = Ti + G with a reciprocal vector G, these generalized Kramers 
pairs are orthogonal to each other: 

((^W(f,)|</^t)(fO)=o, 

since we have 

((0i-*-^(f.)l0i-Hf,)) = e^^"(F')((0(+)(f,)|e,7-'l0i+HfO)) 

= e^^.(F0((e2^-V(+)(f,)|eeHf,))) 



e 



^e„(r,)//^(+) 



(f,)|r,-ie0W(f.))> 



= -e^^'^(^')((eHfoierrVl+Hf.))) 

= -((0(+)(f,O|0t)(f,)). (89) 

We also find that the generalized Kramers partner \<j)n (k)) satisfies (15^ . 

v\4>V{k)) ^ -e^^"(^)e|0W(-fc))> 

= -e'^"(^)ery-i|0l+)(-fc)) 

= -I0^nfe)»- (90) 

Namely, in the presence of the time-reversal symmetry with 0^ = +1 and {rj, 9+} = 0, the states with ((8T|) are 
paired with those with (1821) . As is shown below, this structure enables us to introduce a non-trivial Chern number 
even for the time-reversal invariant system. 

2. Chern numbers for non-Hermittan systems 

Here we generalize the Chern number for non-Hermitian systems. For non-Hermitian systems, the gauge potential 
Ai{k) is introduced as follows [18l.l3q. 

Mk)=Z Y. {{Un{k)\-^\Un{k)), (91) 

ReB„(fc)<B * 

where |u„(fe)) denotes the right eigenstate of H{k) and |w„(fc))) the corresponding left eigenstate. Here we suppose 
that there is a band gap at E in the real part, so no momentum k satisfies Rc£'„(fe) ~ E. The gauge potential Ai{k) 
is well defined in the whole region of the momentum space under this assumption. The Chern number C is defined as 



^-27 



p / / dk^dkyF,y{k), (92) 
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where Fxy{k) is the field strength of the gauge potential Ai(fc), 



^-y^^>- dk, dky ' ^^^^ 

and the integration is performed in the first Brillouin zone. In a manner similar to Hermitian systems, one can show 
that C is quantized and takes only integer values: The Chern number (1921) counts vorticities of wavefunctions in the 
first Brillouin zone [14|. By the gauge transformations \u'^{k)) = e*"^'") |u„(A;)) and \u'^{k)'j} = e*"('^^|u„(fc))), the 
gauge potential Ai{k) is transformed as 

A,(fc)' = A,(fc)-^, (94) 

where the normalization condition (j73p was used. As was shown in Ref.[14], the Chern number C reduces to the line 
integral of the second term of (|94|) around the zeros of the wavefunction, 

^ ^^"(^) dk, (95) 



27r J dk 

which counts the vorticity of the wavefunction in the first Brillouin zone|14|. 

When the system is time-reversal invariant, one can show that the Chern number C is always zero. Thus wc cannot 
use the Chern number itself to characterize topological phases for time-reversal invariant systems. However, if the 
system is pseudo-Hermitian as well and the time-reversal symmetry 0_|_ satisfies Q\ = -t-1 and {77, 64.} = 0, one can 
define a different Chern number which can be non-trivial even in the presence of the time-reversal invariance. The 
key structure is the generalized Kramers pairs explained in the previous section. Under the above assumption for 0+ 
and 77, one can take the following basis, 

r7|</.i±)(fc))=±|<^(±)(fc))), (96) 

with 

(4^)(fc)|0(±)(fc))) = U\^\kMt\k)) = 5^n. (97) 

Here the states \4>n {k)) form the generalized Kramers pair, 

\c^[-\k)) ^ e'^'-C^^e^y-Vi+H-fc))), \A-\k))) = e^'-'^^^Qij-Mt\-k)). (98) 

Using the sign difference of the right hand side of ((96)) . we can introduce the following two different gauge potentials 



A''^\k) = ^ Y. ii'l^tHk)\^\4^\k)}, (99) 

RG£;„(fc)<_E * 

and the corresponding Chern numbers, 

^^^^ = ^ / / dk^dkyFifik), (100) 



/FBZ 

with 

As we mentioned in the previous section, if -E„(fc) is not real, the states |0„ (A;)) are not eigenstates of H{k) in 

general. However, for the real part of the eigenenergy, \(j)n (k)) has the same structure as the eigenstate of H{k). 
Thus the summation in (l99t is well defined. It is also found that 

C = C^+^+C^-\ (102) 

Under the time-reversal transformation, these Chern numbers transform as C^^' — > —C^^K Thus the time-reversal 
invariance implies that C — C'^+-' -|- C*-"^ = 0. On the other hand, the time-reversal invariant combination 

C(+) _ c^-) 
Ctri = ^^^, (103) 

can be non-zero. We call Ctri as time-reversal invariant Chern number. In Sec. llVCSI and Sec. |Vl we will see that 
the time-reversal invariant Chern number characterizes the gapless edge state for a class of non-Hermitian systems. 
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3. Application to 5*0(3, 2) model 



Let us now examine Ctri for the 50(3,2) model (|65|) . Below, we suppose that we have a gap around E^{k) 
and E^{k), and that the system is half-filling. In this model, two right eigenstates |iii (k)) and \u_ (fc)) with the 
negative eigenenergy E_ (k) are 



|uW(fc)) 



a5-E 
^2E{E - as) ( *^i - ^2 



l"^^(fc)) 




^2^(^-05) 



(104) 



where E = \fa\ 



h\ — h\. They form the generalized Kramers pair. 



Here we present explicit calculations of the time-reversal invariant Chern number for 50(3, 2) Hamiltonians ((65 
Substituting ([TUi)) into dMl), we obtain 



Af\k) = 



A\-\k) = 



2E{E - as 



d d 

{E - "'5)t^{E - as) + (as - «a4)^(a3 + «a4) 



(9 



+ {ibi+b2) — {ibi~b2) 



■A 

9fc 



In 



1 



y/2E{E - as) 



2^(^-as) 



9 (9 

(-E - "'5)-^{E - as) + (as + «a4)^(a3 - «a4) 



5 



+ («&i-fe2)^(i&i+62) 



•A 

''dk. 



In 



1 



^2£;(£; - as) 



From ()105p . we find 



i(-) 



i+)i 



A^->ik)^[A^^>i-k)]*, 



(105) 



(106) 



from which we can show C = O'"''-' -|- O'^-' = explicitly. 

To evaluate the time-reversal invariant Chern number, we adiabatically deform the Hamiltonian of the system 
without gap closing in the real part of the bulk energy. This process does not change the time-reversal invariant 
Chern number. In particular, to calculate the time-reversal invariant Chern number for weak non-Hermiticity in 
FigsIH (a) and (b), we decrease the non-Hermiticity (^1,^2) adiabatically as d^ — > (i = 1, 2). In this particular limit, 
we find that the time- reversal invariant Chern number Otri coincides with the spin Chern number Cg in Ref.|ll|. 
Therefore, from the adiabatic continuity, we obtain 



Otri — 2, 



(107) 



for the model in Fig|51 This means that the existence of the gapless edge states in Figl5]is ensured by the non-zero value 
of Otri. Here we should emphasize that the spin Chern number itself is not well-defined once the non-Hermiticity 
(rfij'^2) is turned on. On the other hand, the time-reversal invariant Chern number Otri is well defined even in the 
presence of the non-Hermiticity. 



V. NON-HERMITIAN KANE-MELE MODEL 



So far, we argued systems with imaginary on-site potentials. In this section, we consider a system in which the 
non-Hermiticity is caused by asymmetric hopping integrals. 

Consider the following non-Hermitian version of Kane-Mele model|7l.ll5|: 



H — Hk + Hv + H: 



K -I- nv -I- nso + J^R 



Hr, 



(108) 



where 



Hk = t'^{4cj + H.c), iJy = Ay ^ £,iclci, Hso = i>^so X! (^ii4s''cj -I- H.c), 



i'U) 



ii'U)) 



Hr = -Ar ^ |^c|(s X d^j)^Cj + H.c.j 



(109) 



(^J> 
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FIG. 10: The honeycomb lattice. 
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FIG. 11: The sign of i/i. 



The first term Hk in (|108p is the nearest-neighbor hopping. The second term Hy is sublattice potentials with 
i^i — +1 (—1) for closed (open) circles in FigfTOl The third term i?so represents the next-nearest-neighbor spin-orbit 
interaction preserving the z component Sz of the spin, with vij = ±1. (See Fig llll ) The last term H^, is the imaginary 
Rashba term which gives rise to non-Hermiticity. Here, d = y, di = —^x — iy, and d2 = +^x— iy, are illustrated 
in FigUl 

The imaginary Rashba term Hfi gives asymmetric nearest-neighbor hopping integrals. The term iJ/j is explicitly 
written as 



Hb = -Ar ^ 1^4(5 X dyO^Cj + H.. 






Ar^ [-4-d^xC^ + cl_^U,c, + c\_^U2C^ , (110) 



where J2i, {J2io) denotes the summation over closed (open) circles in FiglTUl Ui = ^<Tx — 2'^v^ ^^'^ ^2 = \<^x + 2'^v 
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FIG. 12: d, di, and da 
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FIG. 13: The nearest-neighbor hopping integrals of imaginary Rashba terms. The hopping integrals of each bond are asym- 
metric, [f/i = \a^ - ^Oy and U2 = \cr^ + ^(Jy] 



The hopping integrals are asymmetric for each bond as illustrated in Fig llSI Let us respectively denote two eigenstates 
of cr^r, Ui, and U2 as |±)^, \±)^, and |±)2: 



a.|±)^ = ±|±)^, C/i|±), =±|±),, t/2|±)2 = ±|±), 



They are explicitly given by 



W"^7!(± 



=j-r 



.ill /Z 



!)• W'^7ir±l J' '^>'^^/2V ±1 



=j-r 



-in jZ 



(111) 



(112) 



For each eigenstate, the hopping integrals t + \rGx-, t + XjiUi, and t + \rU2 read 

(t + Afla.) |±), = (t ± Afl) |±), , (t + A;,C/i(2)) |±)i(2) = (i ± Afl) |±)i(2) , (113) 

and the hopping integrals t — XrOx, t — XrUi, and t — XbU2 yield 

{t - Xna,) |±), = (i T Afl) |±), , (i - XrU,^^)) |±)i(2) = (^ T A^,) \±),^^) . (114) 
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The sign difference between (I113P and (I114|) implies that the hopping integrals (for the basis (|±)^ , |±)-^ , |±)2)) along 
the bonds parallel to (d, di, d2) are asymmetric: (t + A/j) in one direction and (i — A/j) in the other. 

We also note that the imaginary Rashba spin-orbit interaction can be regarded as an imaginary SU{2) gauge 
potential [see Appendix |X] . In this sense, our model can be regarded as a non-Abelian generalization of the Hatano- 
Nelson model in which an imaginary f/(l) vector potential J23l425j was considered. 

Now we consider the Hamiltonian in the momentum space. By performing the Fourier transformation with respect 
to (n, m) in Fig llOi the Hamiltonian in the momentum space is obtained as 

i/(fc)- ^ d„(fc)r,+ Y. C?„6(fc)r,b+Z ^ da(fc)fa+l Y. dab{k)fab, (115) 



a=l,2 a6=12,15 a=3,4 ab=23,24 

where r„b = [Fa, rb]/(2i), and 

di = i(l + 2cosxcosy), d2 — Xv, 

da = Ar(1 — cos X cosy), d^ = — voAr sin a; sin y, 

di2 = — 2icosa;sin?/, (ii5 — Ago (2 sin 2a; — 4 sin a; cos y) , 

(i23 = — ARCOSxsiny, (i24 = VoAr sina;cosy, (116) 

with X = kxa/2 and y — y/3kya/2. Here, we have adopted the following gamma matrices: 

f (1,2,3,4,5) = (o-^^ «) /, o-^ (g) /, cr^ (g) cr^, <T^ (g) cr^, cr^ (g) a^). (117) 

The Hamiltonian (|115p possesses the pseudo-Hermiticity, 

riH{k)^7^-^ = H{k), 7? = zf4f3, (118) 

and time-reversal invariance with Q^ = +1, 

Q+H{-k)Q^^ = H{k), Q+^if^fi-K. (119) 

These symmetries are simply understood by noticing that the Hamiltonian (|115[) is a special case of (|62p by identifying 
(ri,r2,r3,r4,r5) in pTf|) with (F3,F5,F2,Fi,F4) in dnH). From the anticommutation relation {77,6+} = 0, the 
time-reversal invariant Chern number Ctri can be introduced in a manner similar to the previous sections. In the 
following, by using the time-reversal invariant Chern number, we argue the topological stability of the edge states of 
the Hamiltonian (llOSp with zigzag edges [Figl2]. 

First, discuss the Hermitian case Ar = [3, Il3- For Ar = 0, in addition to the time-reversal symmetry 8+ with 
G'^ = +1 (|119p . the Hamiltonian (IllSp has another time-reversal symmetry 0_ with 0^ = —1: 

e-H{-k)ez^ ^ H{k), e-^faFs-X. (120) 

Moreover, for Xn = 0, the pseudo-Hermiticity (|118p reduces to 

[iJ(fc),77]=0. (121) 

Thus the metric operator 77 becomes a conserved quantity whose eigenvalue is either +1 or —1. By regarding rj as 
the z-component of the spin, the spin Chern number Cs can be defined [12] ■ For Ay < 3-\/3Aso, where the spin 
Chern number is non-zero, i.e. Cg — 1, gapless edge modes appear as shown in Fig |14l (a). On the other hand, for 
Ay > S-s/SAso, the system is in the topologically trivial insulating phase with Cs = 0. In this phase, gapless edge 
modes do not appear as illustrated in FiglMlfb). 

Now we include the non- Hermitian term Ar. Once the non-Hermitian term is nonzero, rj is no longer conserved, 
thus the spin Chern number Cg is not well defined. Using the time-reversal invariant Chern number Ctri, however, 
we can argue the topological stability of the gapless edge modes. Let us first consider the region Ay < 3-\/3Aso as 
shown in FigslT5] and [161 When Ar is small, the gapless edge modes that appear in the Hermitian case (Ar = 0) 
still remain [FigslT5](a) and (b)]. These gapless edge states are topologically protected by the time-reversal invariant 
Chern number Ctri — ^'- For small Ar, the non-Hermitian Hamiltonian (|115p can be adiabatically deformed into the 
Hermitian one (Ar = 0) without closing the bulk gap, and the time-reversal invariant Chern number Ctri reduces 
to the spin Chern number Cg in the Hermitian limit. Thus, from the adiabatic continuity, we have Ctri = Cg = 1. 
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A,v=0.1 A.R=0 




FIG. 14: The energy bands of the Kane-Mele model p08|) with zigzag edges along the x-direction for t = 1.0, Aso = 0.06, 
Ar = (Hermitian case), and (a) \v = 0.1, (b)Ay = 0.4. Here a is the lattice constant, and k^ the momentum in the 
a;-direction. In (a), a gapless helical edge state appears on each edge, while in (b) no gapless edge state is obtained. 



On the other hand, for sufficiently large Ar, the bulk gap closes and Ctri becomes ill-defined. Correspondingly, the 
gapless edge modes disappear. 

Next consider the region Ay > 3v3Aso- In the presence of weak non- Hermit icity Ar, gapless edge modes do not 
appear [Figs lTTf a) and irST a')]. This is because the bulk gap has not yet closed and the system remains a topologically 
trivial phase with Ctri = 0. If we further increase Ar, the bulk gap closes near ak^ ~ 27r/3 and Att/3 at a critical 
value of Ar, then the bulk gap opens again [Figs. HD^b) andlWbV] In this region of Ar, gapless edge modes appear 
which are topologically protected by the time- reversal invariant Chern number. For sufficiently large Ar, the bulk 
gap closes once again at ak^ ^ tt, and the edge modes disappear as shown in Figs. [TTT c) and \Wlc) . 



VI. SUMMARY AND DISCUSSION 



In this work, we investigated edge modes and their topological stability in non-Hermitian models. We analyzed 
three types of models: SU{1, 1) lattice model realized on graphene with pure imaginary sublattice potential, 50(3, 2) 
Luttinger model on 2D square lattice, and 5*0(3, 2) Kane-Mele model with asymmetric hopping integrals on graphene. 
The energy spectra of such non-Hermitian models generally contain complex eigenvalues. In this paper, we focused 
on the real parts of the edge bands and characterized them by using topological arguments. (The imaginary part of 
eigenvalues brings decay of wavefunction with time.) For the 5[/(l, 1) lattice model, with numerical calculations, we 
found that the edge states with Re_E = are robust against small non-Hermitian perturbation. We gave topological 
arguments for the robustness of edge state. Meanwhile, the 50(3, 2) Luttinger and the 50(3, 2) Kane-Mele models are 
time-reversal invariant non-Hermitian models with 0^ — +1. The generalized Kramers theorem suggests the existence 
of helical edge modes in the models. The numerical calculations indeed confirmed the existence of helical edge modes 
and robustness of them under the small non-Hermitian perturbations. We introduced time-reversal invariant Chern 
number inherent to non-Hermitian models, and gave topological arguments about the stability of helical edge modes. 

In this paper, we adopted non-Hermitian models whose Hermitian counterparts are typical topological insulators in 
2D. There are many types of topological insulators, such as topological superconductors, 3D topological insulators, etc. 
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FIG. 15: The real part of the energy bands of the non-Hermitian Kane-Mele model (|108p with zigzag edge along x-direction. 
We plot the results for t = 1.0, Aso = 0.06, \v = 0.1, and various values of non-Hermitian parameter Ar. Here a is the lattice 
constant, and kx is the momentum in the a;-direction. In (a) and (b), we have gapless edge states. In (c), the bulk gap closes, 
and no gapless edge state exists. 



It would be interesting to consider non-Hermitian generalizations of these various topological insulators, where gapless 
edge modes could also appear. The topological arguments for non-Hermitian systems presented in this paper would 
constitute the first step to introduce topological invariants characterizing gapless edge modes of other non-Hermitian 
models. 



Note added 



After the completion of this work, we noticed Ref. |37| in which the 5'C/(1, 1) and 5*0(3, 2) models are discussed in 
the context of topological insulators. Their numerical calculations of the lattice version of the SU{\, 1) coincide with 
our results. In Ref. |.37j] . the authors concluded that the appearance of the complex eigenvalues is an indication of the 
non-existence of the topological insulator phase in non-Hermitian models. In the present paper, we focused on the 
real part of the complex eigenvalues and explored the robustness of the gapless edge modes on the basis of topological 
stability arguments. The "topological phase" in the present paper is referred to the phase in which real part of edge 
modes is stable under small non-Hermitian perturbations. 
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FIG. 16: The real part vs imaginary part of the energy bands of non-Hermitian Kane-Mele model (|108p with zigzag edges along 
the x-direction. We plot the results for t = 1.0, Aso = 0.06, \v — 0.1, and various values of the non-Hermitian parameter Ar. 
The data (a) and (b) indicate that the energy of the gapless edge states [Fig llSI (a) and (b)] is real. On the other hand, the 
bulk state supports the imaginary part of the energy. 
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Appendix A: Imaginary Rashba interactions and imaginary 5(7(2) gauge potentials 



Let us first consider continuum Hamiltonians with imaginary Rashba couplings iX (A: real constants), 

H = —{pI + pI) + iXiPxCTy - Pycr^) 

= ■^{Px + iScTy)'^ + —{py-i6(Tr,f + mX'^, 



(Al) 
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FIG. 17: The real part of the energy bands of the non-Hermitian Kane-Mele model pOSp with zigzag edge along s-direction. 
We plot the results for t = 1.0, Aso = 0.06, Xv = 0.4, and various values of non-Hermitian parameter Ar. Here a is the lattice 
constant, and kx is the momentum in the a;-direction. In (a) no gapless edge state exists, while in (b) topologically protected 
gapless edge states appear. In (c), the bulk gap closes and the topologically protected gapless edge states disappear. 



where 9 = mX. From (jAip . we find that the imaginary Rashba interaction can be regarded as an imaginary SU{2) 
gauge potential: 



{A^,Ay) = {-i9ay,i9(7^) . 



(A2) 



These imaginary gauge potential acquires scale transformations Ux and Uy when we proceed by a unit length along 
the X and y directions, respectively: 



Now we perform scale transformations H' = U{x, y) ^HU{x, y) with 
For the first term of the Hamiltonian (jAip , we have 



U{x,y)-^{j)x+i9cryfU{x,y) 



Va^y^-Va^x | _ ■_^ _^ -g^^ | ^Va^x^-Ua^y 



e "e 



(A3) 
(A4) 



„ecr^yJ:^ -ea^y „ 2 
dx^ Px- 



(A5) 
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FIG. 18: The real part vs imaginary part of the energy bands of non-Hermitian Kane-Mele model (|108p with zigzag edges along 
the a;-direction. We plot the results for t = 1.0, Aso = 0.06, Xv ~ 0.4, and various values of the non-Hermitian parameter Ar. 
The data (b) indicate that the energy of the gapless edge states [Fig llTl fbl] is real. On the other hand, the bulk state supports 
the imaginary part of the energy. 



For the second term of the Hamiltonian (|A1I) , we have 

d .„ 



-2ie xya^^-0ayX^0a^y I _j^_ j^Q^ \ ^-9u^y ^Sa-yX ^2i9 xyu^ 



dy 



d^ 



O C CO 

oy^ 

_ _ -2ie-'xytT^ ^ 2ie-'xya^ 

dy-^ 



o \ 2 

— + 2ie'^a,x ] ={py + 29^a^x) ^ , 
where we used the Baker-Campbell-HausdorfF formula up to the leading order of 9. Therefore, we have 



(A6) 



(A7) 
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FIG. 19: Tight-binding model on the square lattice with imaginary 5(7(2) gauge fields. 



The Hamiltonian (|A7|) is Hermitian and coincides with the one with real Rashba coupling A (under the trans- 
formation av — >■ — CTz). We recall that Rashba interactions with real couplings A can be regarded as SU{2) gauge 
potentials [23], 



which give field strengths, 



{A^,Ay) = (-e(jy,ea^). 



i'xy Ox Ay — OyAx — l[AxTAy\ 

= ie'^[(jy,ax\^2e^cT,. 



(A8) 



(A9) 

Actually, (|J7| with a^ ^ -a-, gives field strengths Fxy = 29'^a^ (K9\i . 

Although non-Hermitian Hamiltonian H is transformed into Hermitian Hamiltonian H' by the scale transforma- 
tions U{x,y), non-Hermiticity affects the boundary conditions. Suppose that right eigenfunctions il)'^{x,y) and left 
eigenfunctions tp'^ix, y) of H have periodic boundary conditions both in x and the y directions: 

^«(L„ y) = V^^(0, v), ^^(L., 2/) - ^^(0, y), 

^^{x,Ly)^^^{x,0), 7/;^(x,Lj,) = V^(x,0). (AlO) 

From the boundary condition in the x direction, eigenfunctions tp'i^^ v) — U{x, y)^^ip^{x, y) of H' are found to satisfy 

i^'{Lx,y) - U{Lx,y)-^^\Lx,y) 

= C/(L„y)-i[/(0,y)C/(0,y)-V''(0,y) 
"e-'^--^^'(0,y) 



gScr^y^-eayL.^^-0cr^y^l/( 



— p-ScTyL^ ~2ie a^L^y jl 



V''(0,2/), 



and 



t — ,//('n „^tp2ie <y,L^yJ(7yL^ 



ij'{Lx,yy = ^P'{0,y) 



Similarly, the boundary condition in the y direction gives, 

i^'ix, Ly) = e'^'^^y^'ix, 0), V'(:c, Ly)^ = i^'ix, 0)te-^'^^'^». 

Let us now generalize the above arguments to lattice systems. For simplicity, we consider a tight-binding model on 
the square lattice with the nearest-neighbor hopping integrals with imaginary SU{2) gauge fields: 



(All) 
(A12) 
(A13) 



H = 



Z^ Hn+l,' 



1)6 ^C.{n,m) + C(„_i „)£ 



-'(n^m) 



(n,77i) 



(7i 777 + 11 ^f^77 7T7I ~r L./^, — 1 \ C. C(' 



-(77,m) ~r ''(,i,,77^i)'^ '^(77,771) 

I \^ r t -Sa^ , t ea-y 

+ 1^ [V,m)S C(„+1,™)+C(n,m)e ''C(„_i,„) 



(n,m) 



+*^(n,m)^ "C(n,m+1) + ^{n,m)^ "C(n,m-1) 



(A14) 



which is illustrated in FigfT9l By scale transformations, 
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(n.m) (n.m) ' *• '' 



-{n,m) c e <^{n,m)y ^{n,m) '-■{n,m) 

the Hamiltonian (IA14[) is transformed into 

^ ^ Z^ |_'^(n+l,m)'^(n,™) + '^(n-l,m)'^(n,™) 



/ . |_*^(n,m)*^(^+l,^) "^ *^(n,m)*^(^-l,^) 



(n,7n) 






(A16) 



where we used the Baker-Campbell-Hausdorff formula up to the leading order of 6. The Hamiltonian (JA16P is the tight- 
binding model with SU{2) gauge fields, which gives an observable flux 29^(Jz per plaquette (under a transformation 

Let us impose the periodic boundary condition on the original Hamiltonian H (IAf4|) : 

Then, non-Herniiticity appears as boundary conditions for the Hamiltonian H' (lAf 6[) : 

C{L^,m) — e e C(o,,„), c^^^^^-| — c^p^^^e e , 

C(„,L,) = e^'^^^''C(„,o), cj„,,.^) = cl^of^'^^'-'- (Af8) 
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